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Abstract. - We study many-particle diffusion in 2D colloidal suspensions with full hydrody- 
namic interactions through a novel mesoscopic simulation technique. We focus on the behaviour 
of the effective scaled tracer and collective diffusion coefficients Dt{p)/Do and Dc{p)/ Do, 
where Do is the single-particle diffusion coefficient, as a function of the density of the colloids 
p. At low Schmidt numbers Sc — 0(1), we find that hydrodynamics has essentially no effect on 
the behaviour of Dt{p)/Do. At larger Sc, Dt{p)/Do is enhanced at all densities, although the 
differences compared to the case without hydrodynamics are minor. The collective diffusion 
coefficient, on the other hand, is much more strongly coupled to hydrodynamical conservation 
laws and is distinctly different from the purely dissipative case. 



Introduction ~ The dynamics of Brownian particles in confined geometries, and in two 
dimensions (2D) in particular, is an important theoretical problem with applications in surface 
science and colloidal systems PP . Examples of fundamental questions that have been addressed 
recently are the form of effective interactions between macroions in a colloidal suspension 2 3j 
and the effects of hydrodynamic interactions (His) on the diffusive properties of colloidal 
particles [11310113. 

So far, most studies have focused on self- diffusion of particles in 2D. In the ideal case 
of no external potential and without His, the density dependent self-diffusion coefficients of 
2D hard disk particles have been recently determined using numerical simulations 'H','^ and 
various theoretical approximations O.lOiTTl. While the single-particle limit in the ideal case 
is trivial, no exact analytic results exist for finite densities. In this regime complicated many- 
body effects manifest themselves through memory effects in the motion of tagged colloidal 
particles. 

The situation is even more complicated when the His mediated by the solvent in a colloidal 
suspension are taken into account. Recent work on the self-diffusion of colloidal particles in 
2D and quasi-2D ^USIEIIT] indicate that His do indeed influence self-diffusion. In many cases 
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it has been found that diffusion is enhanced ^|31|H1[7| by the His. Nevertheless, it is possible 
that the nature and magnitude of these subtle effects in a given system depend on the relative 
importance of the His and other interactions. 

While the self-diffusion properties of 2D colloidal systems are relatively well understood, 
much less is known about collective diffusion in 2D colloidal suspensions. Nevertheless, col- 
lective diffusion plays a crucial role in processes such as spreading and phase separation, as 
it describes the decay rate of density fluctuations in a system. While the case without His 
has been considered recently [HlEl i the theoretical understanding of 2D situations with His 
is very limited. This is, in part, due to theoretical difficulties when dealing with collective 
transport in 2D liquids. Furthermore, realistic numerical simulations of collective diffusion 
in two-phase colloidal systems with full hydrodynamics have turned out to be a considerable 
methodological challenge. 

In this letter, our purpose is to employ a recently proposed mesoscopic simulation method 
|13[I14| to shed light on some of the fundamental issues of many-particle diffusion in 2D 
colloids. To this end, we consider an ensemble of colloidal particles that interact mutually 
with a short-range repulsive potential and long-range His. The main issue we want to address 
here is the influence of hydrodynamics in the diffusive dynamics of this system, as the extent 
of His is varied in a controlled fashion by tuning the Schmidt number of the system (see 
below). In particular, we are interested in the behaviour of the effective collective diffusion 
coefRcient as a function of the density of the colloidal particles and the impact of the His on 
this dependence. 

Schmidt number - An important quantity measuring the properties of a fluid in equilibrium 
is the dimensionless Schmidt number iSc, defined as the ratio of momentum diffusivity to mass 
diffusivity: 

Sc^^- (1) 

Here v = rj/ ps is the kinematic viscosity of the fluid, 77 being the viscosity, ps the density of the 
fluid and D the tracer diffusion coefficient of the fluid particles. In a real fluid such as water 
Sc — 0(10'^). Theoretical arguments, too, often include the assumption that hydrodynamic 
fluctuations have reached a steady state on the timescale of the motion of the colloidal par- 
ticles. The situation can be quite different in computer simulations. In molecular dynamics 
(MD) simulations the Schmidt number is typically large. Unfortunately, MD is computation- 
ally too intensive for simulations of colloidal suspensions with explicit solvent, especially at 
high densities. Mesoscopic simulation techniques such as dissipative particle dynamics (DPD) 
reduce the computational cost. On the other hand, in recent DPD simulations Sc^ 0{1) 
It is not quite clear how a Schmidt number several magnitudes smaller than those measured 
in realistic systems influences the properties of e.g. colloidal suspensions. In our opinion, the 
understanding of the influence of Sc on colloidal dynamics is necessary. It is also crucial to 
develop new approaches to bridge the various scales in colloidal systems and, consequently, to 
allow a direct comparison between physical model studies and experiments. 

Model and methods - The model system we consider here comprises an ensemble of disks 
immersed in a 2D liquid. The dynamics of the system is simulated using a novel mesoscopic 
technique introduced by Malevanets and Kapral (MK) |18II14| . The MK method is essentially 
a hybrid MD model, where the colloid is treated microscopically and the solvent obeys coarse- 
grained dynamics. The interaction potential between two solvent particles is zero, while 
the solvent-colloid and colloid-colloid interactions can be chosen at will. For the purposes 
of the coarse-grained solvent dynamics, the time is partitioned into segments r, and the 
system itself is divided into so-called collision volumes or cells. During one collision step t, 
the system obeys Newton's equations of motion, i. e. the time evolution is taken care of by 



E. FALCK et ai: INFLUENCE OF HYDRODYNAMICS ON MANY- PARTICLE ETC. 



3 



Sc 




r [tlj] 


4 M 


Q 


ms [m] 


L[a] 


1 


0.82(1) 


0.5 


2 


±90° 


1 


200 


20 


3.70(2) 


0.1 


2 


±125° 


0.5 


100 


100 


9.11(2) 


0.1 


2 


±170° 


1 


100 



Table 1 - A summary of the three different parameter sets used. See text for details. 



MD. The effective interactions between two solvent molecules occur at each r: a collision 
event takes place, and the velocities of the solvent particles are transformed according to 
Vi{t + t) = V + o; ■ [vi{t) — V]. Here is the velocity of particle i, V is the average velocity 
of all the particles in the cell the particle i belongs to and u; is a random rotation matrix 
over an angle a chosen separately for each cell. It can be shown TS' that this multiparticle 
collision dynamics conserves the momentum and energy in each collision volume. Furthermore, 
it gives a correct description of the hydrodynamics of the velocity field. The method also 
reduces computing times significantly, as the tedious calculation of solvent-solvent forces can 
be omitted. 

The interactions between the colloid particles are strongly repulsive and of short range. 
The interaction potentials are of the form 

I 0, r > Tc, ^ ^ 

where r is the separation between two particles. Here we set n = 12, which allows a direct 
comparison with previous calculations on a smooth surface without hydrodynamics (2!- The 
colloid-colloid interaction parameters were chosen to be CTcc = 2ct and ecc — e, while the colloid- 
solvent parameters are given by CTcs = cr or zero(^) and Ccc = £• The colloidal particles have a 
mass of rric = 5m and the solvent mass is set to = to or to^ = 0.5to. The parameters cr, e 
and to now define our system of units, and hence our unit of time is given by r^j = cr(TO/e)^/^. 

The simulations were performed at a temperature ksT = 2 and for the solvent density 
we used ps = 1. (^) The equations of motion were integrated using the velocity Verlet 
algorithm with a time step At — 0.005. The linear size of our simulation box was set to 
either L = 200 or L — 100. Periodic boundary conditions were employed in all simulations. 
Within our statistical accuracy, the scaled diffusion coefficients did not change with system 
size. However, we did not carry out a systematic finite-size scaling analysis. 

The choice of the parameters that determine the collision dynamics fixes the properties of 
the coarse grained solvent, and in particular the Schmidt number. The value of Sc is affected 
e.g. by the choice of the collision step r, the linear size of the collision volume £c and the 
collision rule, i.e. the choice of the random rotation matrix u). We have used different sets of 
parameters in the simulations. These sets lead to different Schmidt numbers and kinematic 
viscosities. A list of the three different parameter sets is shown in Tabled The systems with 
Sc > 1 have a mean free path (distance traveled during t) of the order of A w O.li'c. In this 

{^)ln the case where the direct solvent-colloid interaction is absent, the colloidal particles participate in the 
collision events llbl . In other words, the solvent-colloid interaction is described indirectly through collisions. As 
a consequence the calculations speed up significantly. We have verified explicitly that the results obtained from 
simulations without solvent-colloid interactions are in quantitative agreement with those computed including 
the interaction directly, see fig.|U(b) 

(^)The solvent and colloidal densities have been defined as ps = Ns / (A/ a"^) and p = N/{A/a'^^), respectively. 
Here A is the area of the system, Ns the number of solvent particles and N the number of colloidal particles. 
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case, to avoid any unphysical correlations at short time and length scales, we have employed 
the random grid-shifting procedure proposed in ref. |17j . 

Results and discussion - Perhaps the simplest transport coefficient is the tracer diffusion 
coefficient Dt, which describes the motion of a tagged tracer particle as Dt = limt^oo Drit), 
where 



among N identical particles. Here d = 2 is the dimension of the system and Vi(t) is the 
velocity of particle i at time t. The quantity (f){t) = (vi(t) •Vi(0)) is the velocity autocorrelation 
function. In evaluating the diffusion coefficients we employ the memory expansion method (18| . 

For colloidal systems with His, it has been known since the 1970s that there are long-time 
tails in the velocity autocorrelation functions. These have been observed in MD J^l and 
lattice-Boltzmann |20| computer simulations, and will lead to a divergence of the tracer diffu- 
sion coefficient in 2D |19I21| . In this case, the diffusion coefficients Drit) can be considered to 
be time-dependent, effective quantities. When tracer diffusion in concentrated suspensions is 
investigated experimentally, it is common to single out the effects of the His by concentrating 
on the so-called short-time diffusion coefficients |22ll4l5ll7| . These are measured at times much 
shorter than t/, the time it takes the tracer particles to diffuse the average distance between 
the suspended colloidal particles. For this definition to be meaningful, such times should be 
significantly larger than the time tb it takes for the velocities of the colloidal particles to 
relax. In our case, for concentrated solutions and small values of Sc in particular, the time in- 
terval between tb and t/ becomes very narrow. Hence, this definition for short-time diffusion 
coefficients cannot be used. 

We consider Dt{p, t) normalised by the single-particle diffusion coefficient Do{t) = Dt{p — > 
0,t). As shown in fig. ^a), in the limit of large times, the quantities DT{p,t)/ Do(t) seem 
to converge to a finite limit. Although the slow logarithmic divergences are present in our 
data (data not shown), the amplitudes of the tails at late simulation times are very small 
and partially masked by statistical fluctuations. In practice we have first determined Dq at a 
time tes where the tracer particle has diffused approximately £ times its own diameter, where 
£ = 25, 5, 2 for Sc — 1, 20, 100, respectively. The effective diffusion coefficients Dt{p) at ffiiite 
densities have then been determined at the corresponding time intervals. 

In fig. n^b) we show a summary of our results for the scaled Dt{p)/Dq as a function of 
the density of the disks p. For comparison we also present our earlier MD results for the 
purely dissipative case Remarkably, it is clear that at low values of Sc hydrodynamics 
has virtually no influence on tracer diffusion. As Sc is tuned to larger values, we find an 
enhancement of Dt{p)/Dq which is largest at intermediate densities. This sheds some light 
on the issues addressed in ref. [S]. The authors study the influence of His on self-diffusion, 
and find that in quasi-2D systems with long-range interactions between the colloidal particles, 
the scaled diffusion coefficients are enhanced. The larger the particle density, the stronger 
the enhancement jSj. Our data indicate that even for short-range interactions of the form 
V{r) ~ 1/r" with n = 12, the His may enhance self-diffusion. Although the mechanism 
leading to enhancement may be complex in a many-particle diffusion process, we expect that 
it is related to the formation of vortices, which have been shown to enhance self-diffusion in 
simple liquids ^ . Further, although direct comparison with experiments is not possible due 
to a different magnitude and range of interactions, we may note that this result is in agreement 
with previous findings ^ElIZI- 

Another important transport coefficient characterizing collective density fluctuations is the 
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Fig. 1 - (a) Scaled effective tracer diffusion coefficients Dripjt) / Do{t) as a function of t for Sc « 1. 
These data have been obtained from simulations without colloid-solvent interactions. The colloidal 
densities from top to bottom are 0.1232, 0.2464 and 0.3697. (b) Scaled effective tracer diffusion 
coefficients Dt{p) / Dq as a function of p for 5c ~ 1 (open circles), Sc ~ 20 (open triangles) and 
Sc K, 100 (open boxes) . For reference, results without His |^ are also presented (solid circles) . 



collective diffusion coefficient Dc{t), which can be defined through a Green-Kubo relation: 

Dc{t) = ^Dut) = e (4) 

Here ^ = {N)/[{N'^) — (N)'^] is the thermodynamic factor (proportional to the inverse of 
the isothermal compressibility kt) and R(i) = ~ is the center of mass 

(CM) displacement '23j. The thermodynamic factor ^, which is a static quantity, is not 
affected by hydrodynamics within our accuracy. To estimate ^ we have used the Boublik 
approximation [22 which is in excellent agreement with MD simulations for most densities [H] . 

In the case of dissipative hard spheres on a smooth surface, it is an exact result that 
the CM mobility Dcm is independent of the density p ^Ol- This is because the interparticle 
interactions preserve the CM momentum, and thus Dcm{p) — ^cm(O) = i^o- However, with 
His in place this argument does not hold any longer. In fig. \^a) we show the scaled CM 
mobility as a function of density in the present system('^). It is now a decreasing function of p, 
and decreases much more rapidly than Dt{p)/Dq shown in fig.^b). It is also remarkable that, 
unlike self-diffusion, Dcmip)/ Dq is not sensitive to Sc: while the self-diffusion of individual 
particles is influenced by Sc, the effects on individual particles are largely independent of each 
other, and thus cancel out in the CM mobility. 

In fig. |2tb) we show the scaled collective diffusion coefficient that displays a minimum at 
small values of density. This minimum is due to the rapid initial decay of Dcm{p)/ Dq. It is 
in striking contrast to the dissipative case (see fig. I3b)) which is entirely determined by ^. 
We argue that this behaviour is generic in colloidal suspensions governed by hydrodynamics, 
since then CM mobility is not constant, but competes with f . However, the actual form of 
Dcip)/ Dq may vary from one system to another, depending on the details of interactions and 
the impact of the His. 

('')Thc CM mobilities Dcm{p) have been determined from Dcm{p,t) in the same manner as the effective 
self-diffusion coefficients Dt{p) have been obtained from DT(p,t). 
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Fig. 2 - (a) Scaled CM mobility Dcm{p) / Dq as a function of p for Sc ~ 1. The inset shows the 
thermodynamic factor ^ from the Boublik approximation ^24j . (b) Scaled effective collective diffusion 
coefficients Dc{p)/ Do as a function of p. Open circles are results with Sc ~ 1, obtained from 
simulations without colloid-solvent interactions. Open diamonds are results from simulations where 
the colloid-solvent interactions are present. Squares, in turn, are from computations where Sc 100. 
For reference, the solid line is Dc{p) / Dq in the case where His have not been taken into account 0. 

Concluding remarks - We have presented a detailed numerical study of the effects of 
hydrodynamics on both self- and collective diffusion of 2D repulsive colloidal particles. This 
has been achieved by using a novel hybrid mesoscopic scheme for two-component liquids. We 
have found that the effective tracer diffusion coefficient Dt{p) / Dq is enhanced by Hfs, an effect 
which becomes more pronounced with an increasing Schmidt number. A more dramatic change 
has been found in the behaviour of the effective collective diffusion coefficient Dc{p)/ Dq, 
which is completely changed by hydrodynamics. This is due to the dynamic term which now 
competes directly with the thermodynamic factor, leading to non-monotonic behaviour in 
contrast with the case of purely dissipative dynamics. It is an interesting question how these 
results are affected when changing from a 2D to a confined 3D geometry as in thin films, and 
how effective interactions obtained from experiments would be manifested in the collective 
diffusion. Work in this direction is in progress. 

* * * 
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